Molecular configurations in a droplet detachment process of a complex liquid 
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We studied the microscopic polymer conformations in the 
droplet detachment process of an elastic semi-dilute poly- 
electrolytic Xanthan solution by measuring the instantaneous 
birefringence. As in earlier studies, we observe the suppres- 
sion of the finite time singularity of the pinch-off' process and 
the occurrence of an elastic filament. Our microscopic mea- 
surements reveal that the relatively stiff Xanthan molecules 
are already significantly pre-stretched to about 90 % of their 
final extension at the moment the filament appears. At later 
stages of the detachment process, we find evidence of a con- 
centration enhancement due to the elongational fiow. 
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I. INTRODUCTION 

The addition of a tiny amount of polymer to a simple 
liquid alters the dynamics of a droplet detachment pro- 
cess dramatically. Instead of the finite time singularity of 
the minimum neck diameter [1], a cylindrical filament is 
formed and the shrinking dynamics can be slowed down 
by several decades [2-5]. Despite the technological rele- 
vance of the droplet forming process of complex liquids 
that reaches from plotting of DNA-microarrays to food 
processing, only little is known on the underlying physi- 
cal mechanisms. Besides the complexity of the problem 
and the lack of a universally applicable constitutive equa- 
tion for complex liquids, there exist no microscopic data 
on the molecular conformation in such a flow. It is the 
goal of this study to close this gap. 

Earlier experimental studies on the droplet detachment 
process of complex liquids were limited to the analysis 
of macroscopic quantities like the determination of the 
shape of the thinning filament. These can be compared 
with theoretical analysis on capillary break-up based on 
phenomenological polymer models like e.g. the Oldryd-B 
or models that follow from kinetic theory like FENE-P 
[6]. But for a true comparison with microscopic models, 
the dynamics of the polymers on the molecular level must 
be measured. Birefringence is a typical tool for such in- 
vestigations [7] , and we use a high speed set-up capable of 
capturing dynamics on the time scale of milliseconds to 
study both the microscopic conformation and the macro- 
scopic fiow response in the elastic thread simultaneously. 



II. SAMPLE CHARACTERIZATION 

The polymeric system of our study, Xanthan (Sigma- 
Aldrich), was chosen because of its high optical 
anisotropy. The molecular weight of t our polymer m 
is vaguely specified by the distributor with "m = 2 x 10^ 
amu or more". Xanthan shows a pronounced shear thin- 
ning and biological digestibility and, for both reasons, it 
is widely used in the food, pharmaceutical, oil, and cos- 
metic industries. The ordered molecule exists in solution 
as a semi rigid helix with a persistence length of 120 nm. 
It undergoes a conformational transition to a disordered 
ficxible coil only above a temperature of 50°C [8,9]. 
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FIG. 1. Stationary viscosities from shear rate and stress 
controlled measurements. Not all experiments and data 
points are shown. The lines are approximations by the Car- 
reau model. 

Our sample solutions were prepared with weight con- 
centrations Cp in the semi-dilute regime from 2bQppm < 
Cp < 2000ppm, mostly dissolved in 80/20 weight% 
glycerol-water as a solvent with a viscosity of rjsoiv = 
dOmPas. We also used different glycerol-water weight 
ratios resulting in solvent viscosities down to rjsoiv — 
2mPas only to allow for a better determination of the 
molecular weight of our sample by fitting our extensional 
rheological data. The polymer concentrations where cho- 
sen to obtain a sufficiently high birefringence but still to 
be significantly below concentrations at which sponta- 
neous lyotropic ordering is expected in equilibrium (cp > 
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5000ppm) [8]. Standard rheological measurements, using 
a cone plate geometry on a standard Rheometer (Haake 
Mars, Thermo Electron, Germany) were performed to 
characterize the shear thinning behavior (fig. 1) of our 
samples. Stationary experiments controlling either shear 
rate or tension were repeated several times for several 
samples in a cone plate geometry with an angle of 0.5° 
and a diameter of 60mm. Depending on the experimental 
deviations an average of between 6 and 10 measurements 
is shown. The solutions show a pronounced shear thin- 
ning, as to be expected for relatively stiff molecules ap- 
proaching a stationary value slightly enhanced compared 
to the pure Newtonian solvent according to the poly- 
mer concentration. The data have been approximated 
using the Carreau model [11] and the results are in rea- 
sonable agreement with ref. [10] if one takes our higher 
solvent viscosities into account. The characteristic relax- 
ation times of the fits range from about 10s for SOOppm 
and 20s for lOOOppm to about 50s for the 2000ppm solu- 
tion. These time scales give a measure for the rotational 
diffusion time of our molecules and must be compared 
with the elongational rates in our droplet experiments. 

For the later discussion of our elongational viscosity 
measurements we would like especially to point out that 
the 1000 ppm samples with different solvent viscosities 
show similar zero shear rate viscosities, and the differ- 
ences only become obvious at high shear rates. 



III. EXPERIMENTAL SETUP 

The experimental setup is sketched in fig. 2A. The 
polymer solutions are quasistatically driven through a 
nozzle of fi = 2mm diameter by a syringe pump. At a 
distance of D = 7cm below the nozzle, a plate is mounted 
in order to omit gravitational effects, and to stabilize 
the filament against air currents. All measurements have 
been performed at room temperature (25°C). The re- 
tardation 6 is mcasiired rising an optical train consisting 
of an 18mW HeNe Laser (JDS Uniphase), a Polarizer P, 
a Photo Elastic Modulator PEM (Hinds Instruments), 
the sample liquid, an Analyzer A, a bright field lens sys- 
tem LI and L2, and a photo diode PD connected to two 
Lock-Ins (Stanford Research). The polarizer is fixed at 
an angle of 45° relative to the photo elastic modulator 
which allows for a Lock-In technique for the detection of 
the retardation signal. The bright field lens system (fig. 
2B) in front of the detector prevents any light that has 
not crossed the filament (fig. 3)) from reaching the detec- 
tor by blocking the parallel components with a mask M 
in the focal point of lens LI. The lens L2 collects the light 
that is diffused by the cylindrical filament into the detec- 
tor [12]. To calculate the birefringence An = 5h{t) / {2'k\) 
(with A = 633nm as the HeNe laser wavelength) from the 
retardation signal the instantaneous thickness h{t) of the 
filament must be known. The thinning process of the fila- 



ment diameter hit) is measured at 500 frames per second 
with a high speed CCD-camera (encore mac PCI lOOOS) 
that is placed perpendicular to the optical train. The 
camera is equipped with a 2 times magnification objec- 
tive and the image capturing is synchronized to the data 
collection of the two Lock-Ins. From fig. 3 it becomes 
obvious that, for geometric reasons, it is impossible to 
measure birefringence before the droplet has completely 
passed the light beam and the temporal exponentially 
thinning filament is present. This is the moment when 
our measurements start. 
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FIG. 2. The experimental setup. See text for further ex- 
planation. 
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FIG. 3. Photographs of the droplet detachment process of a 
2000 ppm Xanthane solution (7730(1, =60mPas) taken every 0.1s 
indicating, schematically, the laser beam positioned close to 
the neck and in a region as cylindrical as possible throughout 
the measurement. 



IV. EXPERIMENTAL RESULTS 

A. Macroscopic measurements 

Different regimes can be distinguished analysing the 
temporal behavior of the neck diameter h{t) in the 
droplet detachment process of a polymer solution. First, 
it follows an exponential law that is well described by the 
Rayleigh instability of Newtonian liquids [13,14]. Then it 
might pass into the self similar regimes that describe the 
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break-up of Newtonian liquids [15] until the polymers in- 
tervene the flow abruptly and a filament is formed which, 
again, thins very slowly and exponentially with time [16]. 
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FIG. 4. Filament diameter h(t) for different polymer con- 
centrations and solvent viscosities. The parameters refer to 
eq. 1. Not all data points are shown. 
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FIG. 5. The apparent elongational viscosity for the same 
set of data as in fig. 4 

Figure 4 shows the filament thickness h{t) extracted 
from the video images that start at the moment when 
the filament is formed. More viscous solvents or higher 
polymer concentrations lead to a slower dynamic that 
differs by more than an order of magnitude. However, 
at first all runs show the same exponential thinning be- 
havior, followed by a linear regime, and the data can be 
approximated by the expression 

h{t) ^ expi-t/r) -bt + yo. (1) 

The formation of the cylindrical filament occurs due to 



the high resistance of the polymers to elongational fiow 
and this resistance is macroscopically quantified by the 
elongational viscosity rj^. The same form of equation 1 
has been used to fit data from extensional rheology mea- 
surements (both capillary break up and filament stretch- 
ing devices [17]). Its physical meaning is that first the 
polymer molecules are uncoiled by the fiow and the dy- 
namic is governed purely by elastic effects. Once the 
polymers are stretched, a steady state is reached which 
leads to the linear dynamic. From h(t) one can calcu- 
late the corresponding elongational rates e = dth{t)/h{t). 
These have to be compared with characteristic relaxation 
times. The microscopic polymer rotational diffusion has 
been estimated by the time constants from the standard 
rheological measurements to be at least A > 10s and the 
appropriate ratio is given by the Weissenberg number 
which is always found to be very large Wi — Xe > 100 in 
our experiments. This means that the polymer molecules 
in the droplet detachment process do not experience a 
force field that is comparable to the Brownian forces 
and time scales, but are affected by a violent flow and 
stretched and oriented in a deterministic manner. Never- 
theless, the elongational rates and the Weissenberg num- 
bers are practically constant in the regime of exponential 
thinning. They diverge, though, when the linear thinning 
behavior sets in. 

The h{t) data can also be used directly to calculate the 
apparent elongational viscosity by equating the capillary 
pressure with the clastic stresses [3,4,17]: 



2a 
th{t) 



(2) 



This approach is a simplification because it neglects 
gravitational or nonlinear elastic effects, but it is com- 
monly used to extract the elongational viscosity from 
capillary break up experiments. A common quantity rep- 
resenting the stress- and stretching history of the fluid, 
and thereby the polymers, is the Henky strain given by 
C = J ^dt. The corresponding elongational viscosity data 
in fig. 5 first show an exponential increase of the ap- 
parent elongational viscosity followed by a steady state 
value for strains C > 1.5. The observation of the plateau 
value at strains that much smaller compared to earlier 
work [4,17,18] is a consequence of the relative stiffness of 
the Xanthan molecule and the pre-stretching before the 
occurrence of the filament. This is also indicated by the 
elevated elongational viscosity at C = 0. If the filament 
would start to be formed with the polymers at rest one 
would expect the elongational viscosity to grow from the 
Newtonian value of the so called Trouton ratio rj^/rj^ ~ 3. 
When rye reaches the plateau value the thinning dynamic 
becomes linear and the elongational rate e diverges, in- 
dicating that is independent of the elongational rate 
as predicted by rigid rod models as well as elastic models 
for Wi » 1 [19]. The measured elongational viscosi- 



February 2, 2008 



4 



ties are in reasonable agreement with data from studies 
with opposing nozzles [20], fiber spinning devices [21], or 
values obtained by analyzing contraction flow [22] 

For our birefringence study we needed an estimate for 
the molecular weight of our sample and we used Batche- 
lors formula [?] for the elongational viscosity of semidi- 
luted rigid rods to fit our data. Fluorescence microscopy 
studies on DNA show that at steady state for Henky 
strains ( > 1.5 the molecules are extended to a maxi- 
mum hydrodynamic length £ that is more than 80% of 
their contour length [23] , and the relatively stiff Xanthan 
molecules are mostly well approximated as rigid rods. In 
order to allow for a robust fitting procedure for the deter- 
mination of the two free parameters molecular weight m 
and polymer length £ we varied both the polymer concen- 
tration and the solvent viscosity. As expected we found 
the elongational viscosity to always be proportional to 
the solvent viscosity (see fig. 5). With our two indepen- 
dent sets of data, we performed the fit with Batchclors 
formula for semidiluted solutions of rigid rods (fig. 6) 

„ ^ hsoiMe/d? 

Tie = orisolv H , 

ln{2£/d) - ln{l + 2./WJ^£/d) - 1.5 (3) 

$ is the particle volume fraction and, for the width of 
the polymer, we took the literature value of d = 2nm. 
The fit yields ^ = A.bfim and m — 2.5 Mamu. This ratio 
of molecular weight and length is in very good agreement 
with e.g. the results from a contraction flow study by 
[22], but differs by a factor of two from values obtained 
by light scattering [20]. The number of Kuhn steps of 
our sample is about 40, a value that is 10 times smaller 
than that of A-DNA which has been used in a previous 
study where no steady state of the elongational viscosity 
had been observed [18]. 
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FIG. 6. The elongational viscosity for different solvent vis- 
cosities (squares) and polymer concentrations (circles). The 
continuous and the dashed lines show the theoretical values 
according to Batchelors formula for semidiluted rigid rods. 

B. Birefringence measurements 

We can now turn to the microscopic study of the con- 
formation of the macromolecules by means of the bire- 
fringence measurements. An elaborate introduction in 
optical flow rheometry is given in [7]. In our unidirec- 
tional stretching geometry with the flow along the di- 
rection of gravity, we can assume transverse isotropy for 
polarizability and segmental order of the single units of 
the macromolecules. Therefore, the birefringence signal 
A7V is proportional to the segmental order of the poly- 
mers and to their concentration 

27r 1 

AN = AaoCpNa- < Scos^ 6* - 1 > (4) 

nm 2 

where n is the average refractive index, Na the Avo- 
gadro number and 6 the orientation angle of the poly- 
mer segments to the direction of extension. By knowl- 
edge of the maximum birefringence signal AN„iax at 
maximum extension and alignment of the polymers, 
the mean fractional extension of our polymers is then 
xjl— ^J[2AN/ ANmax + l)/3 because, on a microscopic 
scale, a variety of different configurations will certainly 
be present [23]. 

The negligence of different configuration types like, 
e.g., dumbbell or hairpin, leads to a certain overestimate 
of the fractional extension and our calculation looses its 
accuracy for coiled states. The correctness of the formula 
improves, though, increasing degree of orientational or- 
der and is a good approximation for large fractional ex- 
tensions of the stretched polymers. 

To test the reliability of our method we performed sev- 
eral runs at different polymer concentrations, and the 
birefringence signal was always roughly proportional to 
the concentration (fig. 7b). The robustness and se- 
lectivity of our method becomes clear if one compares 
the data of the solution with bQQppm of polymer and 
Vsoivent — aOmPas with the solution with lOOOppm of 
polymer and "qsoivent = 15mPas. Both data sets prac- 
tically lead to the same flow (fig. 4) but, exceeding the 
capabilities of simple extensional rheometry, the birefrin- 
gence signal is sensitive to the differences in concentra- 
tion. This shows how important the microscopic mea- 
surements are if one wants to qualify theoretical predic- 
tions that are based on kinetic models. 

The measurements of the birefringence in our setup 
are restricted to the time tc after the filament is formed 
because optical abberations make it impossible to mea- 
sure it at earlier stages of the pinching process. Even if 
the origin of the Hencky strain is not influenced much 
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by this limitation, as most of the strain is accumu- 
lated at the final stages of the experiment, our data 
reveal that the polymers have a significant (flow and 
strain) history when our measurements start. Assum- 
ing that the polymers are completely stretched when 
elastic effects are no longer observable for C > f.5, we 
can calculate the maximum birefringence per polymer 
concentration ANmax/cp = 8 x 10~^g~^cm^ . This is 
comparable to the value for fully elongated A DNA of 
ANjnax/cp ~ 5.1 X I0~'^g~^cm^ [24], but a decade larger 
than values obtained for Xanthan in shear flow [25] and 
in transient electric birefringence measurements [26]. 

The corresponding maximum birefringence of the Cp < 
lOOOppm solutions is only 30% larger than at C = 0. 
This would correspond to a fractional polymer exten- 
sion x/£ = 0.9 at C = 0, indicating that stiff polymers 
like Xanthan are almost completely aligned and stretched 
when the filament occurs. The uncoiling process in the 
range < C < 1-5 correlates with an elastic, exponential 
filament thinning dynamic and a monotonic increase of 
the birefringence signal. 
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FIG. 7. a) The instantaneous birefringence versus the 
elapsed time. The arrows indicate a Hencky strain of = 1.5 
b) the instantaneous birefringence versus the Hencky strain. 
The vertical line indicates the = 1.5 position, the dashed 
lines are linear fits with slopes indicated by the numbers. 
Symbols are as in fig. 1. 

However, instead of a plateau value at the final stage 
of the experiment, we observe a divergence in the bire- 
fringence signal (fig. 7a). The divergence is very repro- 
ducible and takes place at stages of the filament thinning 
process when the filament diameter is still large enough 
to render the analysis unambiguous. The divergence of 
the birefringence signal comes along with a divergence in 
the elongational rate e that follows from the linear shrink- 
ing of the filament. Surprisingly, we find that, the bire- 
fringence signal continues to increase with constant slope 
from C = 1.5toC = 6if plotted against the strain. Only 
for the Cp > IbOOppm solutions did we observe a steeper 
slope in the birefringence signal for < 1.5, presumably 
because of the stronger entanglement of the polymers. 

While the polymers should be completely uncoiled in 
the linear thinning regime, additional physical mecha- 
nisms, that we will discuss in the following, have to be 
considered. We do see three different possible scenar- 
ios: first an effect of the polydispersity of the molecu- 
lar weight distribution of our sample. But the observed 
large Wi should lead to a complete stretching of smaller 
molecules at even smaller strains. Second an overstretch- 
ing of the molecule by the strong flow that would affect 
the configurations of the molecular bonds. Third a con- 
centration enhancement by drainage or evaporation of 
the solvent. This picture is supported by the observation 
that, eventually, the filament might not break at the final 
stages of the thinning process but leave a very thin poly- 
mer (< lO/im) fiber between the nozzle and the ground 
plate. Then a configuration called beads on a string oc- 
curs via an instability of the surface of the cylindrical 
filament [27]. 



V. CONCLUSION 

In conclusion, we have presented the first measure- 
ments on the molecular configurations of semi-rigid Xan- 
than molecules in a droplet detachment process of a com- 
plex liquid. We find that stiff molecules like Xanthan are 
highly oriented and stretched by the flow, even before the 
abrupt transition from the Newtonian self similar shrink- 
ing law to the exponential behavior of the elastic filament. 
Furthermore, we do observe a further increase of the bire- 
fringence signal after the saturation of elastic effects and 
the stretching of the polymers. We discussed different 
possible reasons for this phenomenon and find that a con- 
centration enhancement is most likely to be the pertinent 
effect, but we cannot exclude the possibility that changes 
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in the intramolecular bond configurations caused by the 
violent flow affect the birefringence signal too. 
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